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In systems with a spontaneously broken continuous symmetry, the perturbative loop expansion 
is plagued with infrared divergences due to the coupling between transverse and longitudinal fluc- 
tuations. As a result the longitudinal susceptibility diverges and the self-energy becomes singular 
at low energy. We study the crossover from the high-energy Gaussian regime, where perturbation 
theory remains valid, to the low-energy Goldstone regime characterized by a diverging longitudinal 
susceptibility. We consider both the classical linear O(N) model and interacting bosons at zero 
temperature, using a variety of techniques: perturbation theory, hydrodynamic approach (i.e., for 
bosons, Popov's theory), large- N limit and non-perturbative renormalization group. We emphasize 
the essential role of the Ginzburg momentum scale pa below which the perturbative approach breaks 
down. Even though the action of (non-relativistic) bosons includes a first-order time derivative term, 
we find remarkable similarities in the weak-coupling limit between the classical O(N) model and 
interacting bosons at zero temperature. 
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I. INTRODUCTION 

In the context of critical phenomena, it is well known 
that the Gaussian approximation breaks down in the 
vicinity of a second-order phase transition (below the 
upper critical dimension). When the Ginzburg criterion 
\T — T c \/T c S> t<3 is violated (T c denotes the critical tem- 
perature and \T — T c \/T c ~ to defines the Ginzburg tem- 
perature Tq), the long-distance behavior of the correla- 
tion functions cannot be described by a Gaussian fluc- 
tuation theory and more involved techniques, such as 
the renormalization group, are required (see e.g. (T|). At 
the critical point (T = T c ), one can nevertheless distin- 
guish two regimes in momentum space: a high-energy 
Gaussian regime, where the Gaussian approximation re- 
mains essentially correct, and a low-energy critical regime 
where the correlation function of the order parameter 
field shows a critical behavior characterized by a non-zero 
anomalous dimension r\. These two regimes are separated 
by a characteristic momentum scale pq which defines the 
Ginzburg length £ G = Pq (see e -g- 0)- 

In systems with a broken continuous symmetry, the 
physics remains non-trivial in the whole low-temperature 
phase due to the presence of Goldstone modes, which 
implies that correlations decay algebraically. The cou- 
pling between transverse and longitudinal order parame- 
ter fluctuations leads to a divergence of the longitudinal 
susceptibility |3H5]. Away from the critical regime (i.e. at 
sufficiently low temperatures), one can distinguish a high- 
energy Gaussian regime (|p| 3> Pg)- where the Gaussian 
approximation remains correct, and a low-energy Gold- 
stone regime (|p| -C Pg) dominated by the Goldstone 
modes and characterized by a divergence of the longitu- 
dinal susceptibility. Note that the Ginzburg momentum 
scale pg defined here is the same as the one signaling the 
onset of the critical regime (in momentum space) when 



the system is near the phase transition. For instance, for 
the (if 2 ) 2 theory with O(N) symmetry (classical O(N) 
model), one finds a transverse susceptibility XJ-(p) ~ 
1/p 2 for p — ► 0, while the longitudinal susceptibility 
X\\ (p) ~ l/|p| 4 ~ d is also singular in dimensions 2 < d < 4 
(the divergence is logarithmic for d = 4). At and below 
the lower critical dimension d~ = 2, transverse fluctua- 
tions lead to a suppression of long-range order (Mermin- 
Wagner theorem). There is an analog phenomenon in 
zero-temperature quantum systems with a broken con- 
tinuous symmetry. When the Goldstone mode frequency 
lj = c|p| vanishes linearly with momentum, the longitu- 
dinal susceptibility xy(p,u;) ~ l/(w 2 — c 2 p 2 )( 3 ~ d ^ 2 has 
no pole-like structure but a branch-cut for d < 3, and 
the dynamical structure factor exhibits a critical contin- 
uum above the usual delta peak S(u> — c|p|) due to the 
Goldstone mode (SHE]- 

Historically, the divergence of the longitudinal sus- 
ceptibility was encountered (although not recognized as 
such) early on in interacting boson systems. The first 
attempts to improve the Bogoliubov theory of superflu- 
idity [9] were made difficult by a singular perturbation 
theory plagued by infrared divergences [TU1 - H3"] . As real- 
ized later on [14ITT5] . the singular perturbation theory is 
a direct consequence of the coupling between transverse 
and longitudinal fluctuations. 

In this paper, we study the crossover from the high- 
energy Gaussian regime to the low-energy Goldstone 
regime in the ordered phase, both for the classical 
O(N) model and interacting bosons at zero tempera- 
ture. Even though the action of (non-relativistic) bosons 
includes a first-order time derivative term, which pre- 
vents a straightforward description in terms of a classical 
0(2) model, we find remarkable similarities in the weak- 
coupling limit between these two models. On the other 
hand, the strong-coupling limit of the O(N) model, i.e. 



2 



the critical regime near the phase transition, has no direct 
analog in zero-temperature interacting boson systems. 

The classical 0(A) model is studied in Sec. [TTJ while 
superfluid systems are discussed in Sec. |III| First, we 
show that the loop expansion about the mean-field so- 
lution is plagued with infrared divergences and deduce a 
perturbative estimate of the Ginzburg momentum scale 
(Sees. II A and III A). Then, we use symmetry ar- 
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guments to derive t 



re exact value of the self-energies at 

and 



vanishing momentum (and frequency) (Sees. II A 3 
III A 3 1 . In the case of bosons, we obtain Nepomnyashchii 
and Nepomnyashchii 's result about the vanishing of the 
anomalous self-energy [T4"] . 



II B and IIIB 



we 



In Sees 

show that the difficulties of perturbation theory can be 
circumvented within a hydrodynamic approach (i.e., for 
bosons, Popov's theory |17H19j ) based on an amplitude- 
direction representation of the order parameter field. 
This yields the correlation functions in the hydrody- 
namic regime defined by a characteristic momentum scale 
Pc ^> Pg- The O(N) model is solved in the large- N limit 
in Sec 
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This allows us to obtain the longitudinal 
correlation function in the whole low-temperature phase, 
including the critical regime in the vicinity of the phase 
transition. Finally, we show how the non-perturbative 
renormalization group (NPRG) provides a natural frame- 
work to understand the ordered phase of the O(N) model 
and the superfluid phase of interacting bosons (Sees. II D 
and|ffl~Cl. 



II. THE {ifi 2 ) 2 THEORY AT LOW 
TEMPERATURES 



we have introduced the "normal" (£ n ) and "anomalous" 
(San) self-energies. In the following, we assume that the 
order parameter cp Q is along the direction (1, 0, • • • , 0) so 
that 



E tl (p) 



£ n (p) + £ an (p) if * = 1, 
E n (p) - E an (p) if i ^ 1. 



(5) 



The anomalous self-energy £ an is related to the spon- 
taneously broken 0(A) symmetry and vanishes in the 
high-temperature phase. E n and £ an are analogous to 
the normal and anomalous self-energies which are usu- 
ally introduced in the theory of superfluidity |201 I21j . 
For N — 2, we can introduce the complex field 



1 
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[y>i(r) + iip 2 {r)]. 



(6) 



Making use of the two-component field 



*(r) = 



*t(r) = (V*(r),V(r)), (7) 



the two-point propagator becomes a 2 x 2 matrix in 
Fourier space, whose inverse is given by 



p 2 +r + £ n (p) 
E an (p) 



E an (p) 
* 2 1 - 1 £„(p) 



ro 



(8) 



and bears some similarities with the single-particle prop- 
agator in a superfluid (Sec. Ill I. 



We consider the (<p ) theory defined by the action 

%>] = /^j^v^ + ^ + l^ 2 ) 2 } (1) 

where (p is a TV-component real field and d the space 
dimension. We assume N > 2 and d > 2. The model 
is regularized by a ultraviolet momentum cutoff A. The 
connected propagator 

G iS (p) = (^(pV,(- P )} - (Vi(p)X^(-p)) (2) 

is related to the self-energy E by Dyson's equation G _1 = 
Gq 1 + E, where 



G o^(p)^ p2+ro 



(3) 



is the bare propagator. In the low-temperature phase, 
if we denote by cp = (v( r )) the order parameter, the 
self-energy 

Ejj(p) = £ 0ii ¥> 0)j E ; (p) + (6i,j - £ 0ji £ 0ji )£ t (p) 

= kA^n(p) - E an (p)] + 2y. 0i ^ 0)J .E an (p) (4) 

(<p a = ip /\ip \) can be written in terms of its longitudinal 
(£;) and transverse (E t ) parts. In the second line of Q, 



A. Gaussian approximation and breakdown of 
perturbation theory 

Let us begin with a dimensional analysis of the action 
(fTl. If we assign the scaling dimension 1 to momenta 
(I7e. [p] = 1), the field has engineering dimension [cp] = 
[ro] = 2 and [uq] = 4 — d. We can then define two 
characteristic length scales, 



-1/2 



l/(d-4) 



(9) 



In the critical regime of the low-temperature phase (£ 3> 
Cg)j £,g is the characteristic length scale associated to 
the onset of critical fluctuations, while £ = £j is the 
Josephson length separating the critical regime from a 
regime dominated by Goldstone modes |22| . When criti- 
cal fluctuations are taken into account, one finds that £j 
diverges with a critical exponent v which differs from the 
mean-field value 1/2. At low temperatures away from the 
critical regime (£ <C £g), £ = £c corresponds to a correla- 
tion length for the gapped amplitude fluctuations while 
direction fluctuations are gapless due to Goldstone's the- 
orem. The physical meaning of the Ginzburg length 
in this temperature range will become clear below. 
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FIG. 1. One-loop correction E' 1 ' to the self-energy. The 
dots represent the bare interaction, the zigzag lines the order 
parameter <po, and the solid lines the connected propagator 
G<°>. 



1. Gaussian approximation 

Within the mean-field (or saddle-point) approxima- 
tion, one finds tp — \ip \ — (— 6ro/ito) 1//2 m the low- 
temperature phase (rrj < 0). In the Gaussian approx- 
imation, one expands the action to quadratic order in 
the fluctuations cp — ip This yields the (zero-loop) 
self-energy 



V<WfWi - J _3r ° if i = 1 > 
Hi (P) - i _ rQ if i ^ lf 



(10) 



for |p| <C A, where 



A d = 



2 1 - 



- d / 2 T(d/2) 



1 



sin(Trd/2) T(d-\) 



if d < 4, 
if d = L 



(14) 



The one- loop correction (12 1 diverges for p — > and the 



perturbation expansion about the Gaussian approxima- 
tion breaks down. By comparing the one-loop correc- 
tion to the zero-loop result, i.e. |En (p)| ~ E n °^(p) or 



I San (p) | ~ sin' 1 (p) , one can nevertheless extract a char- 
acteristic (Ginzburg) momentum scale, 

f [A d (N -l) Uo ] 1/(4 - d) if d<4, 
PG ^iAexp( 3Iw ^) if d = 4, < 15 > 

which was obtained previously from dimensional analysis 
[Eq. (|9|]. While the Gaussian or perturbative approach 
remains valid for |p| pc, the limit |p| <C pg cannot be 
studied perturbatively. We shall see in Sec. |II B| that the 
breakdown of perturbation theory is due to the coupling 
between transverse and longitudinal fluctuations. 



from which we obtain the longitudinal and transverse 
propagators, 



Gf (p) = Gff(p) = 



1 

P 2 + 2|r |' 



Gi 0) (p) = G^( P ) = ^. 



(11) 



In agreement with Goldstone's theorem, the transverse 
propagator is gapless, whereas the longitudinal suscepti- 
bility Gi(p = 0) = 1/| 2r*o | is finite. We shall see below 
that this last property is an artifact of the Gaussian ap- 
proximation. 



2. One-loop correction and the Ginzburg momentum scale 

The one-loop correction E^ to the self-energy is shown 
diagrammatically in Fig. [T] While the first diagram is fi- 
nite, the second one gives a diverging contribution to En 
in the infrared limit p — > when d < 4. The divergence 
arises when both internal lines correspond to transverse 
fluctuations, which is possible only for En. Thus E22 is 
finite at the one-loop level and the normal and anomalous 
self-energies exhibit the same divergence, 



EWfpl^Wfp),-^ 



q q 2 (p + q) : 



r, (12) 



where we use the notation J q = J ' momen tum 



integration in ( 12 ) gives [23] 

f 1 = / A d \p\ d - 4 if d < 4, 
lq 2 (p + q) 2 I Ailn(A/|p|) if d = 4, 



(13) 



3. Exact results for E n (p = 0) and E an (p = 0) 

Although the one-loop correction E^'(p) diverges 
when p — > for d < 4, it is nevertheless possible to 
obtain the exact value of E(p = 0) using the O(N) sym- 
metry of the model. 

Let us consider the effective action 



r[0] = -lnZ[h] + Jd d rh-4> 



(16) 



defined as the Legendre transform of the free energy 
— In Z [h] where h is an external field which couples lin- 
early to the ip field and 



5hi(r) 



(17) 



The notation (• • -)h means that the average value is com- 
puted in the presence of the external field h. satisfies 
the equation of state 



= hi(r) 



(18) 



At equilibrium and in the absence of external field, the 
order parameter ip = (<p(r)) is obtained from the sta- 
tionary condition of the effective action, 



ST[d>} 
8ct>i(v) 



= 0. 



(19) 



<M r )=¥>o 



T[(j>] is the generating functional of the one-particle irre- 
ducible vertices 



r (n) fn ••■ r 1 



niri) ■ • • 6<j) in (r n ) 



(20) 



The later fully determine the correlation functions. In 
particular, the two-point vertex is related to the 
propagator by = G^ 1 = G^ 1 + E. 

The 0(N) invariance of the action ([TJ implies that the 
effective action T[<p] is invariant under a rotation of the 
field 4>. Let us consider the case N = 3 for simplicity (the 
following results are easily extended to arbitrary N) . For 
an infinitesimal rotation — > <fi + On x about the axis 
n (n 2 = 1 and 9 0), the invariance of the effective 
action yields 



ijk 



(21) 



where eyfc is the totally antisymmetric tensor. Taking 
the first-order functional derivative 8/5(f)i(r') and setting 
4>i(v) — 5j, t x<po, we obtain 



With n = (0, 0, 1), this gives 



r( 2 2 ) (p = 0)=r + S 22 (p = 0)-0, 



(22) 



(23) 



where r( 2 )(p) denotes r^ 2 )(p,-p). Equation ([23] 
is a direct consequence of Goldstone's theorem. If 
we now take the sec ond -order functional derivative 
<5( 2 )/<^(r')<^ m (r") of plj and set &(r) = fc.i^o, we 



obtain the Ward identity 

E[ r ^( r '. r ")%7+rf (r",r')e ijm 

+ I ddr E r ^( r < r '' r ")^i W = 0. (24) 
Choosing I = 2, m = 1 and j = 3, this gives 

rgV,r") - if^rV) -<p Q j d d rTfl{v,T',v") = 0. 

(25) 

Integrating over r' and r" and using ( 23 1 , we deduce (in 
Fourier space) 



rg>(o,o,o) = r " (0 '°> 



(26) 



where V is the volume of the system. 

Let us now consider the exact diagrammatic represen- 
tation of the self-energy shown in Fig. [2] We know from 
perturbation theory that the third diagram in Fig. [2] is 
potentially dangerous when the two internal lines corre- 
spond to transverse fluctuations. We therefore write the 
self-energy Sn(p) as 



Eii(p) = En(p) 



N - 1 



6VV 



u Q fo ^2 G 22(q)G22(p + q) 



x Iw 2 (-p ) -q,p + q), 



(27) 





r( 4 ) 




FIG. 2. Exact diagrammatic representation of the self-energy 
in terms of the three- and four-leg vertices and I^ 4 '. The 
dots represent the bare interaction, the zigzag lines the order 
parameter, and the solid lines the (exact) connected propaga- 
tor. 



where En(p) denotes the part of the self-energy which 
is regular in perturbation theory (i.e. the part that does 
not contain pairs of lines corresponding to G 22 G 22 ). If 
we assume that the transverse propagator G 22 (q) is pro- 
portional to 1/q 2 for q — > (this result will be shown 
in the following sections), the integral JqG^q) 2 is in- 
frared divergent for d < 4. To obtain a finite self-energy 
Eii(p = 0), one must require that 



Iimrg ) 2 (0 J -q,q) 

q— >0 



r( 3 2 ) 2 (o,o,o) = o. (28) 



The Ward identity (26| then implies T$(p = 0) = 0, so 
that we finally obtain 

£„(p = 0) = -r + \ [r n (p = 0) + r 22 (p = 0)] = -r , 



E an ( P = 0) = - [r n (p = o) - r 22 ( P = 0)] = 0. 



(29) 



It may appear surprising that the anomalous self-energy, 
which is related to the spontaneously broken O(N) sym- 
metry, vanishes for p — 0. The equivalent property in 
interacting boson systems is a fundamental result of the 



theory of superfluidity (Sec. Ill I 



B. Amplitude-direction representation 



The difficulties of the perturbation theory of Sec. |II A| 
can be avoided if one uses the "good" hydrodynamic vari- 
ables in the low-temperature phase, namely the ampli- 
tude and the direction of the ip field. We thus write 



ip(r) = p(r)n(r), 
where n(r) 2 = 1, and obtain the action 



(30) 



S[p,n] = J d<v{l(V P ) 2 + ^(Vn 



2 , r o 2 , "o 4 



4! 



(31) 
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At the mean-field level, the amplitude takes the value 
po = (— Gro/uo) 1 / 2 in the low-temperature phase (ro < 
0) . For small amplitude fluctuations p' = p — p (which is 
expected to be the case at sufficiently low temperatures) , 
we obtain the action 



S\p',n] 



d d r 



\ro\P 



Pi 



(Vn) 2 (32) 



and deduce that the amplitude fluctuations are gapped, 



(p'(p)p'(-p)) 



1 



Pc 



(33) 



If we are interested only in momenta |p| <C p c — \/^\ r a\- 
to first approximation we can ignore the higher-order 



terms in p' that were neglected in (32 1, since they would 



only lead to a finite renormalization of the coefficients of 
the action S[p',n] [23"] . 

Equation (32) shows that in the "hydrodynamic" 
regime |p| <C p c direction fluctuations are described by 
a non-linear sigma model. It is convenient to use the 
standard parametrization n = {a, tt) where a is the com- 
ponent of n along the direction of order and tt a (N — 1)- 
component field (n 2 = a 2 + tt 2 = 1). Integrating over a, 
one obtains 



S[ P ',Tt} 



d d r -W) 



'\2 



\ro\P 



12 



pI 



(Vn) 2 (34) 



for small transverse fluctuations tt [24 . In this limit, we 
can treat TTjfr) as a variable varying between — oo and 
oo. From p4|), we deduce the propagator of the it field, 



(7Ti(P>i(-p)> = 



2 J 

PoP 2 ' 



(35) 



Again we note that the terms neglected in (34) would 
only lead to a finite renormalization of the (bare) stiff- 
ness Pq of the non-linear sigma model at sufficiently low 
temperature. In fact, equation (34 1 gives an exact de- 
scription of the low-energy behavior |p| <C p c if one re- 
places Pq by the exact stiffness and p~ x — (2|ro|) -1 / 2 by 
the exact correlation length of the p' field. 

We are now in a position to compute the longitudinal 
and transverse propagators using 

<pt = pa = p\fl - tt 2 ~ po + p' - \po^ 2 , / 3g x 

ft = P™ — Pan- 

Since the long-distance physics is governed by transverse 
fluctuations, we have retained in (36 1 the leading contri- 



butions in tt. Making use of (35 1, one readily obtains 

1 



G t (p) ^ Po(^(PK(-P)) = 
The longitudinal propagator is given by 

G i (r) = (p'(r)p'(0)) + ^ 2 ( 7 r(r) 2 7r(0) 2 ) c 



(37) 



(P'(r)p'(0)) 



N - 1 

~2P§ 



: G t (r) 



(38) 



where (•• -) c stands for the connected part of (•■•). The 
second line is obtained using Wick's theorem. In Fourier 
space, this gives 



1 



Pc 



N - 1 



1 



q 2 (p + q) 



(39) 



where the momentum integral is given by ( 13 1 for |p| <C A 
and d < 4. By comparing the two terms in the rhs 



of (39 1, we recover the Ginzburg momentum scale (151. 
ForJp| 3> pa, the longitudinal propagator Gi(p)~ 
l/(p 2 + p 2 ) is dominated by amplitude fluctuations and 
we reproduce the result of the Gaussian approximation. 
On the other hand, for |p| < p G , <3j(p) - l/|p| 4_d 
is dominated by direction fluctuations and diverges for 
p^ 0. 

The divergence of the longitudinal propagator is a di- 
rect consequence of the coupling between longitudinal 
and transverse fluctuations [3 . In the long-distance limit, 
amplitude fluctuations become frozen so that \ip\ = p ~ 
Pq. This implies that the longitudinal and transverse 
components tpi and ip t cannot be considered indepen- 



dently as in the Gaussian approximation (Sec. II A) but 
satisfy the constraint <pf + tpf ~ p\. To leading order, 



i Pl ~ p (l - \f/ 2 and G,(r) ~ G t (v) 2 [Eq. (38))], i.e. 
G/(p) ~ 1 / 1 p> | 4 rf for d < 4 (the divergence is logarith- 
mic for d = 4). 

Equations p7| and ( 39 1 imply that the self-energies 
must satisfy 



Sii(p) = -r -p 2 + C 1 |p| 4 - 
£22 (p) = -r + C 2 p 2 , 



(40) 



for p — > and d < 4, i.e 
£ n (p) = -r 

S an (p) 

For d = 4, one finds 

S n (p) = -ro 

San(p) 



+ o(p 2 ), 



^|p| 4 - d + 0( P 2 ). 



(41) 



Ci 
ln(A/|p|) 



0(p 2 



(42) 



ln(A/|p|) 



+ 0(p 2 ). 



For p = 0, we reproduce the exact results of Sec. |II A3| 
Equations ( 41|42 1 show that S n (p) and E an (p) contain 
non-analytic terms that are dominant for p — > 0. 



C. Large- iV limit 

In this section, we show that the previous results for 
the longitudinal propagator are fully consistent with the 
large- N limit of the (y 2 ) 2 theory. Furthermore, the large- 
TV limit enables to compute the longitudinal propaga- 
tor not only at low temperatures but also in the critical 
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regime near the transition to the high-temperature (dis- 
ordered) phase. 

To obtain a meaningful large- N limit, we write the 
coefficient of the (y 2 ) 2 term in Eq. (JlJ as uq/N and take 
the limit N — > oo with uq fixed. Following Ref. [55], we 
express the partition function as 

V[<p, P , A] e -fd d r[H^) 2 + r -iP+i^P 2 + ^^ 2 -P)}_ 

(43) 

It can be easily verified that by integrating out A and 
then p, one recovers the original action S[ip\. If, instead, 
one first integrates out p, one obtains 

Z = J T>[ip, A] e-J^Mi^^il^+^/Afa-ro)^ 

(44) 

As in Sec. |IIB| it is convenient to split the ip field into a 
a field and a (N — l)-component field n. The integration 
over the tt field gives 



(det<?) 



(iV-l)/2 



where 



g-\r, r') = -V 2 5(r - r') + iA(r)5(r - r') 



(45) 



(46) 



is the inverse propagator of the TTi field in the fluctuating 
A field. We thus obtain the action 



S[a,X] = i J d d r [{Vaf+iXa 2 } 



3N 
2u 



d d r {iX - r ) 2 + 



N -1 



Trlng- 1 . (47) 



In the limit N — > oo, the action becomes proportional to 
N (this is easily seen by rescaling the a field, a — > %/iVcr) 
and the saddle-point approximation becomes exact. For 
uniform fields <r(r) = a and A(r) = A, the action is given 

by 



~S[a,X] = ^° 3 



^(iA-ro^ + ^Trlns- 1 (48) 

ZUq lv 



(we use N — 1 ~ N for large A), with <7 _1 (p) = p 2 + iX 
in Fourier space. From ( |48| , we deduce the saddle-point 
equations 



crm = 0, 



6iV 



(m 2 — r ) — iV 



(49) 



where we use the notation m 2 = iX (iX is real at the 
saddle point) . These equations show that the component 
a of the ip field which was singled out plays the role of 
an order parameter. 

In the low-temperature phase, a is non-zero and m = 
0. The propagator <?(p) = 1/p 2 is gapless, thus identi- 
fying the TTi fields as the N — 1 Goldstone modes associ- 
ated to the spontaneously broken O(iV) symmetry. From 



Eq. ( 49 ) , we deduce 



6iV 
u 



where 



roc = 



uo 
' 6 



(ro - roc), 



u K d 



6 d-2 



A 



d-2 



(50) 



(51) 



(with K d = 2 1 - d ir- d / 2 /r(d/2)) is the critical value of 
ro. Since the saddle-point approximation is exact in the 
large- N limit, the effective action T[a, X] is simply given 
by the action S[a,X] [Eq. p7|] gS]. We deduce 



r (2) (l 



r( 2 j(r-r') r^(r-r') 
ri 2 J(r^r') ri 2 A } (r-r') 



(2), 



-V 2 <5(r-r') 
icrS(r - r') f 



iaS(r — r') 
n(r-r') + ^(r-r') 



where 



n(r - r') = g{r - r')g(r' - r) 



(52) 



(53) 



and we use the notation r^o^r — r') = S^T /Sa(r)Sa(r'), 
etc. The two-point vertex is computed for the 

saddle-point values of the fields a and A. In Fourier space, 
we obtain 



r (2) ( P ) 



p 

ia 



IG 

w -n( P ) . 



2 ' «o 

and the propagator G = r' 2 ' -1 takes the form 



G(p) 



1 



detr( 2 )(p) 



f n( P ) + 



with 



detr^(p) = p 2 



N 



n( P ) 



3N 



3N 
u 



-ia 

„2 



(54) 



(55) 



(56) 



and II(p) = / g(q)g(p + q). Equation (561, together 
with the small p behavior of II(p) [Eq. ( 13 FfTleads us to 



introduce three characteristic momentum scales, 



PG = 



PJ 



u A d 



2a 2 



l/(4-d) 



l/(d-2) 



NA d/ 

2\ 1/2 



12 

u A ( 



(roc - r ) 



l/(d-2) 



, (57) 



37V 



[2(r 0c -r )} 1/2 . 



For simplicity, we discuss only the case d < 4; equiva- 
lent results for d = 4 are easily deduced. The Josephson 
length £,j = Pj 1 - which separates the critical regime 
from the Golstone regime (see below) [55] - diverges at 
the transition with the critical exponent v = X/(d— 2). 
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FIG. 3. Characteristic momentum scales pc p.j and p c vs 
T c - T for fixed u [Eqs. (57| with r = Fo(T - T )]. 
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FIG. 4. Momentum dependence of the longitudinal correla- 
tion function G Iir (p) in the critical and non-critical regimes 
of the low-temperature phase as obtained from the large- iV 
limit (2 < d < 4). 



which also characterizes the divergence of the correlation 
length in the high-temperature phase |23| . The momen- 
tum scales ( 57 ) are not independent since 



Pc 



Pi 



— ) 

PgJ 



d-2 



(58) 



If we vary ?'n with uq fixed, we find that the three char- 
acteristic scales (57 1 are equal when T = Tq, where Tq 
is defined by 



ro(T c -T G ) 



upA 
6 



2/(4-d) 



(59) 



(see Fig. [3]). We have assumed that rrj = fo(T — To) 
(with T the mean-field transition temperature) and used 
roc = ^o(T c — To). We recognize in (59 1 the Ginzburg cri- 
terion |2] so that we can identify Tq with the Ginzburg 
temperature separating the critical regime near the tran- 
sition from the non-critical regime at sufficiently low tem- 
peratures. 

In the critical regime (T c - T < T c - T G or pj < 
Pg), using pj <C p c <C pg, °ne finds the longitudinal 
correlation function 



P 2 ,- d 



Gcrcrip) = 



if 



A _ d „ |P |«PJ, 

if IpI >PJ, 



(60) 



while in the non-critical regime (T c — Tc; -C T c — T or 

PG <Pc), 

C(p)=j p ' }l p ^ (61) 

— j if |p| >_p G . 

P 2 +p 2 

In the non-critical regime, we recover the results of sec- 
tion III Bl We find two characteristic momentum scales 
(pc and p c ) and two regimes for the behavior of G aa (p) : 

i) a Goldstone regime (|p| <C pc) characterized by a 
diverging longitudinal propagator Go-o^p) ~ l/|p| 4_rf , 

ii) a Gaussian (perturbative) regime (|p| ^S> Pa) where 
G CTCT (p) ~ l/(p 2 + Pc)- The critical regime is char- 
acterized by two momentum scales (jpj and po) and 
three regimes for the behavior of G aa (p): i) a Gold- 
stone regime (|p| <C pj) with a diverging longitudinal 
propagator, ii) a critical regime (pj <C |p| <C pc) where 
Gacr(p) ~ 1/ |p»| 2 17 with a vanishing anomalous dimen- 
sion r\ (ij is C(l/A^) in the large- TV limit |2.3I 126]). iii) 
a Gaussian regime (p G <C |p|) where G CTCr (p) ~ 1/p 2 - 
These results are summarized in figure [4] 



D. The non-perturbative RG 

1. The average effective action 

The strategy of the NPRG is to build a family of the- 
ories indexed by a momentum scale k such that fluctu- 
ations are smoothly taken into account as k is lowered 
from the microscopic scale A down to |27l 125] . This is 
achieved by adding to the action Q the infrared regula- 
tor 



p.' 



(pi(-p)R k (p)cpi(p). 



(62) 



The average effective action 

T k [</>] = - In Z k [J] + J d d rJ2 Ji & ~ [4>] (63) 

i 

is defined as a modified Legendre transform of — In Z k [J] 
which includes the subtraction of AS k [4>\- Here Ji is an 
external source which couples linearly to the <pi field and 
<p(r) — ((f(r))j. The cutoff function R k is chosen such 
that at the microscopic scale A it suppresses all fluctua- 
tions, so that the mean-field approximation Ta[^)] = S[<f>] 
becomes exact. The effective action of the original model 
([!} is given by T k=0 provided that R k =o vanishes. For a 
generic value of k, the cutoff function R k (p) suppresses 
fluctuations with momentum |p| < k but leaves unaf- 
fected those with |p| > k. The variation of the average 
effective action with k is governed by Wetterich's equa- 
tion VM 



d t T k [cj>] = ^Tr 



Rk ( 



(64) 
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(2) 

where t — ln(fc/A) and R k — dtR k - F k '[<p] denotes the 
second-order functional derivative of Ffc [(f)]. In Fourier 
space, the trace involves a sum over momenta as well as 
the internal index of the cf> field. 

Because of the regulator term AS k , the vertices 



(Pi 



, p„) are smooth functions of momenta 



and can be expanded in powers of pf/k 2 . Thus if we 
are interested only in the long distance physics, we can 
use a derivative expansion of the average effective ac- 
tion |27l I28| . In the following, we consider the ansatz 



r fc [0] = J d d r^{V4>y 



Mp) 



(65) 



Because of the O(iV) symmetry, the effective potential 
U k (p) must be a function of the O(N) invariant p = <p 2 /2. 
To further simplify the analysis, we expand U k (p) about 
its minimum /Oo.fc, 

U k (p) = U k (p Q . k ) + y(p - Po.kf- (66) 

We consider only the ordered phase where po.fc > 0. 
In a broken symmetry state with order parameter <p = 
(^/2 j oo,fc, 0, • • • ,0), the two-point vertex is given by 



r&P) 



Z k p 2 + 2X k p , k if i = 1, 



if i ^ 1. 



(67) 



(2) 

By inverting T k , we obtain the longitudinal and trans- 
verse parts of the propagator, 



G k ,i(p) 



1 



Z k p 2 + 2X k po, k ' 



(68) 



Since these expressions are obtained from a derivative 
expansion of the average effective action, they are valid 
only in the limit |p| <C k. In practice however, one can 
retrieve the momentum dependence of G k= o(p) at finite 
p by stopping the RG flow at k ~ |p|, i.e. G k= o(p) ~ 
Gfc^| p |(p), where Gfc^,| p |(p) can be approximated by the 
result of the derivative expansion. It is possible to ob- 
tain the full momentum dependence of the correlation 
functions in a more rigorous and precise way, within the 
so-called Blaizot-Mendez-Weschbor scheme [30-32], but 
this requires a much more involved numerical analysis of 
the RG equations. 

The transverse propagator Gfc jt (p) is gapless [Eq. ( 68 1], 



in agreement with Goldstone's theorem, which is a mere 
consequence of the O(N) symmetry of the average effec- 
tive action (65 1. On the other hand, the divergence of 



the longitudinal susceptibility obtained in the previous 
sections suggests that X k — > for k — > (limfc_j.o Po,k > 
in the ordered phase). We shall see that this is indeed 
the result obtained from the RG equations. 



2. RG flows 

It is convenient to work with the dimensionless vari- 
ables 

Po,k = Z k k 2 ~ d p ^ k , 



X k — Z, 2 k d 4 Afc. 



(69) 



The flow equations f or pp, k , X k and Z k are obtained by 
inserting the ansatz ( 65|66 l into the RG equation (64 1. 
The calculation is standard [27, 28] and we only quote 
the final result. 



d t po,k 
d t ~Xk 



(2- 
(d 



Vk)Po,k 
2r] k )X k 



3- 
2 h ' 1 



N 



—Ik,t, 



-Xl{9J k , n (0) + (N-l)J k>tt (0)}, 
i lk =2Xl~p 0tk [Jl u (0) + J'u ti (0)], 



(70) 



where rj k = — dtlnZ k denotes the running anomalous 
dimension. With the cutoff function R k (p) = Z k (p 2 — 
fc 2 )9(p 2 -fc 2 ) (33] {Q(x) is the step function), the thresh- 
old functions appearing in ( 70 1 can be calculated analyt- 
ically (see Appendix [A| . 

In Fig. 5] we show A&, i] k and p Q k vs — t = ln(A/fc) for 
d = 3 andiV = 3. We fix Afc = o = wo/3 and vary ro (i.e. 
p 0i fc = o — — 3ro/ito)- When the system is in the ordered 
phase away from the critical regime (red solid lines in 
Fig. [SJ, i.e. p c 3> Pgi we see a crossover for k ~ pq 
(to = ln(pc/A) ~ —4) from the Gaussian regime to the 
Goldstone regime characterized by A& ~ A*, rj k = and 
Po,k ~ k" 1 (i.e. p , k ~ Pq = lim fc ^ Po,fc)- Since A fc ~ A* 
and rj k = imply Afc ~ k, we find that the longitudinal 
susceptibility Gfc : ;(p) = l/2X k p _ k ~ 1/k diverges when 
k — > 0. Identifying k with |p| to extract the momentum 
dependence (as explained above), we recover the singular 
behavior Gfc = o,;(p) ~ VIpI m three dimensions. More 
generally, for an arbitrary dimension, one finds Afc ~ k e X* 
and Gfc,/(p) - l/k e = l/\p\ e with e = 4 - d. Thus 
in the RG approach, the divergence of the longitudinal 
susceptibility is a consequence of the existence of a fixed 
point for the dimensionless coupling constant Afc. 

When the system is in the critical regime of the ordered 
phase (blue dotted lines in Fig. [5|, i.e. pq pj, there 
is a first crossover from the Gaussian regime to the crit- 
ical regime for k ~ pq followed by a second crossover 
to the Goldstone regime for k ~ pj. In the critical 



regime p G > k > pj, X k 



A* 



Vk 



?f and p . k ~ p* 



are nearly equal to their values at the critical point be- 
tween the ordered and disordered phases 1341 [551 . This 
gives Gfc jt (p) ~ G fci /(p) 



l/fc-""p 2 , i.e. Gfc =0 ,t(p) 



Gfc=o.;(p) ~ 1/IpI v ^ we identify k with |p| 



3. Analytical solution in the low-temperature phase 

In the low-temperature phase (away from the critical 
regime, i.e. when p c 3> Pg)- it is possible to obtain an 
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FIG. 5. (Color online) Xk, rjk and po.fc vs —t = ln(A/fc) for 
d = 3, JV = 3, A = 1 and A fe=0 = 0.5. The (red) solid line 
corresponds to po,fc=o = 0.3 (p c 3> pa) and the (red) dots are 
obtained from the analytic solution (731. The (green) dashed 
line corresponds to po,k=o = 0.072147 (p c ~ pc ~ pj) and 
the (blue) dotted one to p ,fc=o = 0.072146123 (pc > pj). 



find 



eA c p c 



tKpl 



ek' 



l^(N-l)X c pl 



(72) 



for k <C p c - The last expression can be rewritten in the 
more insightful form 



A fc = 



A* 



l + ik/paY' 



where 



and 



lim Ai. = ; r 

k^o 8v d {N - 1) 



(73) 



(74) 



Pg 



(N-l) 
(N-l) 



dc 

8vdXc 
deZl 

Pc 



1A 



(75) 



Equation ( 73 ) is in remarkable agreement with the nu- 
merical solution of the flow equations ( 70 ) (Fig. |5| . In the 
weak-coupling limit pq <C p c , we can ignore the renor- 
malization of Z k as well as that of A& between k = A and 
k = p c , and approximate Z Pc ~ 1 and A c ~ A^ = a = uq/3. 
We then recover the expression 



Pg 



(N-l) 



8v d u 



3de 



(76) 



of the Ginzburg momentum scale obtained in previous 
sections. A similar analysis can be made for the case 
d = 4. 



analytical solution of the flow equations for k <C p c - In 
this limit, the RG flow is dominated by the Goldstone 
modes and the contribution of the longitudinal mode can 
be omitted. This amounts to ignoring Jk,u(0), J' k lt (0) 



and J' k t ,(0) in Eqs. (70 I, which is justified by noting that 

Afc/5o,fc becomes very large for k <C p c (X^kPo,k ~ k 2 ~ d for 
k 0), where the hydrodynamic scale p c is defined by 
2A Pe po,p c ~ 1- This gives rjk — and 



dt.X 



(71) 



where v d = [2 d+1 n d / 2 T(d/2)}- 1 . We have used the ex- 
pression of the threshold functions given in Appendix |A"| 
Equation (711 should be solved with the boundary con- 
dition Afc = A c for k — Aq ~ p c . For d < 4, we then 



III. INTERACTING BOSONS 

We consider interacting bosons at zero temperature 
with the (Euclidean) action 



S = I dx 



(77) 



where ip(x) is a bosonic (complex) field, x — (r, r), and 
J dx = Jq dr J d d r. r € [0, /3] is an imaginary time, /3 — > 
oo the inverse temperature, and fj, denotes the chemical 
potential. The interaction is assumed to be local in space 
and the model is regularized by a momentum cutoff A. 
We consider a space dimension d > 1. 
Introducing the two-component field 



*(p) = 



ip(p) 
ip*{-p) 



&(p) = (r(p),i>(-p)) (78) 
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(with p = (p, iuj) and u> a Matsubara frequency), the one- 
particle (connected) propagator becomes a 2 x 2 matrix 
whose inverse in Fourier space is given by 

iu + fj, - e p - E n (p) -E an (p) \ 

~ s an(p) ~ iuj + M - e P - S„(-_p) y ' 

where E n and E an are the normal and anomalous self- 
energies, respectively, and e p = p 2 /2m. If we choose 
the order parameter (ip(x)} — y / no to be real (with no 
the condensate density), then the anomalous self-energy 
E an (p) is real. 

To make contact with the classical (if 2 ) 2 theory with 
O(N) symmetry studied in Sec. [Tl] it is convenient to 
write the boson field 



= 7/| W' 1 0*0 +ii J 2(x)] 



(80) 



in terms of two real fields ipi and ip2 and consider the 
(connected) propagator Gij{x,x') — (ipi(x)ipj(x')} c . The 
inverse propagator G" 1 ^) reads 



(81) 



where 



(82) 



e p - /i + En (p) u + E12 (p) 
-cj + E 2 i(p) e p - [i + E 22 (p) I ' 



E u (p) = ~[E n (p) + E„(-p)] + E an (p), 
E 22 (p) = ~[E n (p) + E n (-p)] - E an (p), 
E ia (p) = |[S n (p) - £„(-?)], 
£21 (p) = --[E n (p) - £ n (-p)], 
when E an (p) is real. 

A. Perturbation theory and infrared divergences 

1. Bogoliubov's theory 

The Bogoliubov approximation is a Gaussian fluctu- 
ation theory about the saddle point solution ip(x) = 
y/™o = \/ H*/ g (i-e. ipi(x) = \J1n§ and ^>i(x) = 0). It is 
equivalent to a zero- loop calculation of the self-energies, 

^(p) = 2gn , S£)(p)= 5 n 0) (83) 

or, equivalently, 

Eg ) (p) = 3 9 n 0l Eg ) (p)=ffn , (84) 
This yields the (connected) propagators 

G<?>(p) = -^(?#*(p)>c = 
G<£?(p) = -Mi#(-p)> £ 



w 2 
gn 



(85) 



where _E p = [e p (e p + 2gno)] 1//2 is the Bogoliubov quasi- 
particle excitation energy. When |p| is larger than the 
healing momentum p c = (2gmn ) 1 ^ 2 , the spectrum E p ~ 
e p + grig is particle-like, whereas it becomes sound-like 
for |p| « p c = \/2mc with a velocity c = yjgno/m. 
In the weak-coupling limit, no — n ( n is the mean 
boson density) and p c can equivalently be defined as 
p c = (2gmn) 1 ' 2 . In the hydrodynamic regime |p| <C p c , 



(0) _ 2gn 

U 22 IP J — 2~~i 2~ 2' 

O 5 ) = 2~~i 2^2 • 



(86) 



Note that in the Bogoliubov approximation, the occur- 
rence of a linear spectrum at low energy (which implies 
superfluidity according to Landau's criterion), is due to 
E an (0) being nonzero. 



2. Infrared divergences and the Ginzburg scale 

Let us now consider the one-loop correction £W to the 
Bogoliubov result E^°\ For d < 3, the second diagram of 
Fig. [T] gives a divergent contribution when the two inter- 
nal lines correspond to transverse fluctuations, which is 
possible only for En. Thus E 22 is finite at the one-loop 
level and the normal and anomalous self-energies exhibit 
the same divergence, 

£«(p) ~ £«(p) ~ ~ l -g 2 n I G^(q)G^{p + q), (87) 



where we use the notation q — (q, iuj') and J — 
f°° f . For small p, the main contribution to the 

J— oo 2-7T Jq 

integral in ( |87[ ) comes from momenta |q| < p c and fre- 
quencies |u/| < cp c , so that we can use (86 1 and obtain 



E« (?)-£«(?) 



4. 3 

,5 no 



Q 



Q 2 (Q 



where Q = (q, uj'/c) and P = (p,w/c) are (d + 1)- 
dimensional vectors. The momentum integral in ( |88| 
< and is given by (13), with A 
1 and |p| by (p 2 + lu 2 /c 2 )^ 2 . 



is restricted by |Q| ^ p c 
replaced by p ci d by d 



We can estimate the characteristic (Ginzburg) momen 
turn scale pq below which the Bogoliubov approximation 
breaks down from the condition |En (p)\ 

\p) \ ~ Si2(p) for |p| = p G and |w| = cp G , 



|E«/ 



si 0) (p) or 



(Aj+^m^/O-d) if d < 3, 
p c exp | - ) if d= 3. 



£ 2 ' 



This result can be rewritten as 

Pc (Ad+ig d/2 ) in3 - d} if d < 3, 



PG 



p c exp 



(89) 



(90) 
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where 



g = gmn 



l-2/d 



\ni/d) 



(91) 



is the dimensionless coupling constant obtained by com- 
paring the mean interaction energy per particle gn to the 
typical kinetic energy 1 jmf 2 where f ~ fT 1 ' 4 is the mean 
distance between particles |36| . A superfluid is weakly 
correlated if g <C 1, i.e. pc *C p c *C n x l d (the charac- 
teristic momentum scale n 1 ^ does however not play any 
role in the weak-coupling limit) [37 . In this case, the 
Bogoliubov theory applies to a large part of the spec- 
trum where the dispersion is linear (i.e. |p| < p c ) and 
breaks down only at very small momenta |p| < pq <C 
p c . In the next sections, we shall see that the weakly- 
correlated superfluid bears many similarities with the or- 
dered phase of the classical O(N) model away from the 
critical regime. When the dimensionless coupling g be- 
comes of order unity, the three characteristic momen- 
tum scales pc ~ p c ~ n 1 ^ become of the same order. 
The momentum range [pc,Pc\ where the linear spectrum 
can be described by the Bogoliubov theory is then sup- 
pressed. We expect the strong-coupling regime g ^> 1 to 
be governed by a single characteristic momentum scale, 
namely n 1 ^. 



3. Vanishing of the anomalous self-energy 

The exact values of E n (p = 0) and E an (p = 0) can be 
obtained using the U(l) symmetry of the action, i.e. the 
invariance under the field transformation ip(x) — > e^ipix) 
and ip* (x) — > e~ i8 ip* (x) [HUSH]. The derivation is similar 
to that of Sec. Ill A 31 Let us consider the effective action 



= - \nZ[J x , J 2 ] + J dxlhfa + J 2 <h], 



(92) 



where Jj is an external source which couples linearly to 
the boson field ipi, and <pi(x) = (ipi(x))j the superfluid 
order parameter. The U(l) symmetry of the action im- 
plies that r[0] is invariant under a uniform rotation of the 
vector field (<j>i(x), <?i> 2 (x)) T . For an infinitesimal rotation 
angle 9, this yields 



If we now take the second-order functional derivative 
5™ /8<pi(y)S<f> m (z) of ( 93 1 and set <pi(x) — 8 i} iy/2no, we 



obtain the Ward identity 

i i 

-y/2no / dxT^ m {x,y, z) = 0. 



(96) 



Integrating over y and z and setting I = 2 and m = 1, 
we deduce (in Fourier space) 



lf 22 (0,0,0) 



1 1^(0,0) 



(97) 



where we have used ( 95 ) 



The self-energy En can be written as 

9\j ^7 ^2 ^2 2 (g)G 22 (p + q) 



x r i22(-p 5 -q,p + q), 



(98) 



where En(p) denotes the regular part of the self-energy 
(i.e. the part that does not contain pairs of lines corre- 
sponding to G 22 G 22 ). If we assume that the transverse 
propagator G 22 (g) ~ l/(w 2 + c 2 q 2 ) at low energies (this 
result will be shown in the following sections), the inte- 
gral J G 22 (g) 2 is infrared divergent for d < 3. To obtain 
a finite self-energy En(p = 0), one must require that 

r i22 (o,o,o) 
ri?(p = o) = 



= 0. The Ward identity (97 1 then implies 
and in turn 



E n (p = 0) = 

E an (p = 0) = 



rg>(p = 0)+r£(p = 0) 



(2), 



r^(p = o)-r$( P = o) 



= o. 



(99) 



The vanishing of the anomalous self-energy E an (p = 
0) was first proven by Nepomnyashchii and Nepom- 
nyashchii [T3] . To reconcile this result with the existence 
of a sound mode with linear dispersion, the self-energies 
E n (p) and E an (p) must necessarily contain non-analytic 
terms in the limit p — >• (Sec. Ill B 4 1. 



dx 



€ij4>j(x) = 0, 



(93) 



where e,j is the totally antisymmetric tensor. Taking 
the functional derivative S/6<f>i(y) and setting (f>i(x) — 
Si i\/2no leads to 



r( 2) ( P = o) = o. 



(94) 



For 1 — 2, equation (94 1 yields the Hugenholtz-Pines the- 
orem [T2] 

lf 2 } (p = 0) = E n (p = 0) - E an (p = 0) - n = 0. (95) 



B. Hydrodynamic approach 

It was realized by Popov that the phase-density rep- 
resentation of the boson field -0 = y/ne t9 leads to a the- 
ory free of infrared divergences |16fFTH] . Popov's theory 
bears some similarities with the analysis of the (if 2 ) 2 
theory based on the amplitude-direction representation 
(Sec. II B |. In this section, we show how the phase-density 
representation can be used to obtain the infrared behav- 
ior of the propagators G n (p) and G an (p) without encoun- 
tering infrared divergences |19] , Our approach is simi- 
lar to that of Popov (with some technical differences in 
Sec. |IIIB2|. 
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1. Perturbative approach 

In terms of the density and phase fields, the action 
reads 



S[n, 9} 



dx 



in9 + (V60 2 

2m 



(Vn) 2 

8mn 



9 2 
2 



(100) 

At the saddle-point level, n(x) — n — u/g. Expanding 
the action to second order in 5n = n — fi, 9 and V#, we 
obtain 

(Vn) 2 



S[Sn, 



dx 



i5n9 



2m 



(V9Y 



fj 



8mn 



[Sri) 



(101) 

The higher-order terms can be taken into account within 
perturbation theory and only lead to finite corrections of 
the coefficients of the hydrodynamic action ( 101[ ) [18] . 

We deduce the correlation functions of the hydrody- 
namic variables, 

G nn [p) = (Sn[p)Sn[-p)) = 



G n6 [p) = (Sn[p)6[-p)) 
Gee{p) = (9[p)6[-p)) = 



J 2 + El ■ 



(102) 



+ a 



whe re E v is the Bogoliubov excitation energy defined in 
Sec. IIIA1 In the hydrodynamic regime |p| < p c = 
y/2gmn, 



G nn (p) 



m uj 2 + c 2 p 2 



G n e[p) = — o" 

0J Z 



c 2 p 2 



(103) 



Gee{p) 



mc 



n lo z + c z p 



2 n 2 : 



where c = gn/m is the Bogoliubov sound mode veloc- 
ity [p c = V2mc). 



where x' = (r',r'). n(x), Vn(x) and id T 9+ ^(W) 2 are 
Galilean invariant (but d T n(x) is not). V 2 6* is also in- 
variant but is odd under time-reversal symmetry. Thus, 
to second order in derivatives, the most general effec- 
tive action compatible with Galilean invariance and time- 
reversal symmetry reads 



T\n, 



,1,4 ^M(Vn) 2 



8m 



U[ri) 



2 

P =i 



Cp(n) 



i<9 T 6> 



2m 



(V6 



(105) 



up to an additive (field-independent) term. K(n), U[ri) 
and Cp(n) are arbitrary functions of n. 

To determine c p [n), we now consider the system in the 
presence of a fictitious vector potential [Aq, A), 



S[, 



;4J = 



in(0 T - A ) + 7— (V6» - A) 2 
2m 



(Vn) 2 
8mn 



/Lin 



f 2 

2" 



(106) 



The action is invariant under the local U(l) transfor- 



mation 9 —> 9 + a and —> A 



n -r d^a where a[x) 



is an arbitrary phase. By requiring that r[n, 0;A M ] 
r[n, 6 + a; A M + d^a] shares the same invariance, we de- 
duce 



r[n,B;A M ] 



d.rA 5^(Vn) 2 + C/(n) 



2 

■^Cp(n) 



8m 
id r — iAq 



2m 



Noting that 



n(x) 



S\nZ[J n , Jg; A^] 
iSAo(x) 



[V9-Af 



ST[n,9;A^] 
iSAq[x) 



(107) 



(108) 



we must have Ci(n) = n and c p (n) = for p > 2. We 
conclude that 



2. Exact hydrodynamic description 



In this section, we show that equations ( 103 1 are exact 



in the low-energy limit |p|, \uj\/c <C p c provided that c 
is the exact sound mode velocity and h the actual mean 
density (which may differ from fj./g). Let us consider the 
effective action T[n, 9] defined as the Legendre transform 
of the free energy — lnZ[J„, Jg] [J n and Jg are external 
sources linearly coupled to n and 9) [3S] . At zero temper- 
ature, T[n, 8] inherits Galilean invariance from the action 
(100 1. In a Galilean transformation (in imaginary time), 



r + ivt and t' — t . the fields transform as 
n[x) = n[x), 



9'[x') = 9{x) 



-rnv t — mv • r, 



(104) 



Tin, 



,/.,,<; ^M(Vn) 2 



8m 



id T t 



\-U[n) 

(Vfl) 2 1 

2m 



(109) 



to second order in derivatives. 



From ( 109 1, we obtain the two-point vertex in constant 
fields n[x) = h and 9[x) — const (with n the actual boson 
density) , 



r$(p) T%[ P ) 



T {2) (P) = I (2), v F (2), 

= ( Y -Bv 2 + U"{n) 



(110) 
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By inverting T^(p), we recover the propagators 
(103 1 in the small momentum limit |p| <C p c = 
[4mC/"(n)/y(n)] 1 / 2 but with a sound mode velocity c 
given by 



nU"(n) 



m 



(111) 



Noting that the compressibility k — n 2 dn/du can also 
be expressed as [30] 



1 



n 2 U"{nY 



(112) 



we conclude that the Bogoliubov sound mode velocity c 
is equal to the macroscopic sound velocity (mn/t) -1 ' 2 . 
Moreover, since the superfluid density n s is defined by 
r 22 } (P:°) = ^P 2 for p ->■ 0, we find that at zero 
temperature n s = n is given by the fluid density [13J. 



3. Normal and anomalous propagators 

To compute the propagator of the ip field, we write 

1>{x) = ^n + Sn(x)e ie(x) , (113) 



where uq = \(-ip(x))\ 2 = \(^n(x)e l9 ^ xS, }\ 2 is the conden- 
sate density. For a weakly interacting superfluid, no ~ n, 
and we expect the fluctuations Sn to be small. Let us as- 
sume that the superfluid order parameter (tjj(x)) = ^/no 
is real. Transverse and longitudinal fluctuations are then 
expressed as 



H2 
fyi 



Sn 
/2n~o~ 



n _0n2 



(114) 



where the ellipses stand for subleading contributions to 
the low-energy behavior of the correlation functions. For 
the transverse propagator, we obtain 



G 22 (p) ^ 2n Gee{p) 



2nomc 2 



n uj 2 + c 2 p 2 
to leading order in the hydrodynamic regime, while 

UJ 



(115) 



2 J_ ^2p2 



Gi 2 (p) ~ G ne (p) 
The longitudinal propagator is given by 
Gn(aO 



(116) 



^- o G nn (x) + ^-(e(x) 2 6(0) 2 ) 



1 

2n 



G nn (x) + n G e e(x) 2 , 



(117) 



where the second line is obtained using Wick's theorem 
(which is justified since the Goldstone (phase) mode is 



effectively non-interacting in the hydrodynamic limit). 
In Fourier space, 



G n (p) 



2mno uj 2 + c 2 p 2 



+ n Gee*G e e(p), (118) 



where 



Gee-kG ge (p) = / Gee(q)G e e{p + q) 



(119) 



with the dominant contribution to the integral coming 
from momenta |q| < p c and frequencies \oj'\/c < p c . Us- 
ing ( 13 1, we find 



Gee * Gge(p) 

A d+1 c(f) 2 ( P 2 + £) 



(d-3)/2 



if d < 3, 
if d = 3. 



(120) 



By comparing the two terms in the rhs of ([118]) with 
|p| = pq and \w\ = cpa, we recover the Ginzburg scale 
(891. For |p|,|a;|/c pc, the last term in the rhs of 
(118 1 can be neglected and we reproduce the result of 
the Bogoliubov theory (noting that h ~ no), while for 

| P |, M/C « p G , G n (p) ~ 1/{UJ 2 + C 2 p 2)(3-d)/2 ig domi _ 

nated by phase fluctuations. The longitudinal suscepti- 
bility Gn(p,ioj = 0) ~ l/|p| 3 ~ d for p —> in contrast to 
the Bogoliubov approximation Gn(p, iuj = 0) = l/2mc 2 . 

From these results, we deduce the hydrodynamic be- 
havior of the normal propagator, 



G n (p) 



~[Gr 11 (p)-2iG ia (p) + G 22 (p)] 

n^mc 2 1 
h u) 
iuj I 



c 2 p 2 



OJ' 



c 2 p 2 



Gu(p), 



(121) 



as well as that of the anomalous propagator, 
1 



G an (p) 



[Gii(p) - G 22 (p)} 



n^mc 



n or + crp 



(122) 



where Gn(p) is given by (118). The leading order terms 
111 ( fl2l| ) and ( fl22| ) agree with the results of Gavoret and 
Nozieres [T3] and are exact (see next section). The con- 
tribution of the diverging longitudinal correlation func- 
tion was first identified by Nepomnyashchii and Nepom- 
nyashchii [T5], and later in Refs. [HH |4"TH4"4"] . 
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4- Normal and anomalous self -energies 

To compute the self-energies £ n (p) and S an (p), we use 
the relations 



E n (p) =G„ 1 (p) 



Gn(-P) 



Sanb) 



G n (p)G n (- P ) - GUpY 
G an (p) 



(123) 



G n (p)G n {~p)-G an ( P y 



with 



G n {p)G n {-p) - G an (pf 

= G 11 (p)G 2 2(p) + G 12 ( P ) 2 



G 22 (p) 



n Gg e *Ge9(p) + 



2riQmc 2 



Setting 



G n (p) ~ --G 22 (p), 
G an (p) ~ ].G 22 {p), 



(124) 



(125) 



in the numerator of Eqs. ( 123 I, we obtain 

EanbHEn^-G^b) 



r (w 2 + c 2 p 2 ) (3 ^ )/2 if d<3, 



111 



( <?pI \ 

\ui 2 +c 2 p 2 J 



if d = 3, 
(126) 



in the infrared limit |p|,|oj|/c <C PGj where G 1 (p) 



e p +/i. Equations (126) agree with the exact results 



(99 1 and show that E n (p) and £ an (p) are dominated by 
non-analytic terms for p — > 0. This non-analyticity re- 
flects the singular behavior of the longitudinal correlation 
function 



Guip) 



1 



2Sa„(p) 



(127) 



in the low-energy limit. 

It should be noted that the singularity of the self- 
energies is crucial to reconcile the existence of a sound 
mode with a linear dispersion and the vanishing of the 
anomalous self-energy £ an (p = 0) |14| . In the low-energy 
limit. 



S an (p) = AS(p) + S an (p), 
Z n (p)-Go\p) = A£(p) + £ n (p), 



(128) 



where A£(p) denotes the singular part (126) while £ n (p) 
and E an (p) are regular contributions of order p 2 ,w 2 . Us- 
ing AS(p) » £ n (p) - Go l (p),^ an (p) for p -> 0, by in- 
verting (|79[) we obtain 



G„(p) 
G an (p) 



1 



2[£ n (p)-£ an (p)]' 
1 

2[E n (p)-£ an (p)]' 



(129) 



Since both E n (p) and £ an (p) can be ex panded to order 
p 2 , uj 2 , we conclude that equations ( |129 1 predict the exis- 
tence of a sound mode with linear dispersion. Of course, 
Eqs. (1291 are nothing but our previous equations (1151 
and (Il25f . 



In deriving the low-energy expression ( 126 1 of the self- 
energies, we have assumed that the hydrodynamic de- 
scription holds up to the momentum scale p c and ignored 
the contribution of the non-hydrodynamic modes. In 
Popov's original approach |19| . one introduces a momen- 
tum cutoff po satisfying pq <g; Pq <C p c - Since p Q 3> pq, 
modes with momenta |p| > pq can be taken into account 
within standard perturbation theory (see Sec. Ill A). On 
the other hand, low-momentum modes |p| < po <C p c 
are naturally treated in the hydrodynamic approach dis- 
cussed in this section. The final results are independent 
of po. The only difference with our results (1261 is that 
p c in the expression of E an (p) for d = 3 is replaced by a 
smaller momentum scale 1451. 



C. The non-perturbative RG 

The NPRG approach to zero-temperature interacting 
bosons has been discussed in detail in Refs. |8l l43l 144} I4BT 
I51| . Our aim in this section is to briefly summarize the 
main results |52| while emphasizing the common points 



with the classical 0(A) model studied in Sec. II D 



To implement the NPRG, we add to the action an in- 
frared regulator term 



(130) 



which suppresses fluctuations with momen- 
tum/frequency below a characteristic scale fc but 
leaves high momentum/frequency modes unaffected. 
The average effective action is defined as 

r fc [^,0] = - In Z k [J*, J] + Y^[J*(p)<t>(p)+c-c-} 

p 

-AS k [ct>*,4>}, (131) 

where 4>(x) — (ip(x))j is the superfluid order parameter. 
J denotes a complex external source that couples linearly 
to the boson field. satisfies the RG equation (64 1. 



As in Sec. |IID| we choose the cutoff function R k such 
that all fluctuations are suppressed for fc = A (so that 
Ta[0*,0] = S[(f>*,(j>}) and Rk=o(p) = 0. In practice, we 
take [S] 



Rk{p) 



2m 



fc 2 



1,2 „2 



(132) 



where r(Y) = (e Y — 1)~ 1 . The fc-dependent variable 
k is defined below. A natural choice for the veloc- 
ity Co would be the actual (fc-dependent) velocity of the 
Goldstone mode. In the weak coupling limit, however, 
the Goldstone mode velocity renormalizes only weakly 
and is well approximated by the fc-independent value 
c = y/gn/m. 
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1. Derivative expansion and infrared behavior 

The infrared regulator ensures that the vertices are 
regular functions of p for |p| -C k and \ui\/c <C k even 
when they become singular functions of (p, icu) at k — 
(c = Cfc ~ Cfc = o is the velocity of the Goldstone mode). 
In the low-energy limit |p|, \uj\/c <C k, we can therefore 
use a derivative expansion of the average effective action. 
We consider the ansatz 



r k [<f>*,(/>] = / dx 



Zc,kd T 
A fe , 



V Ak d 2 



^A,k, 

2m 



+ -^-( n - n o,kY 



(133) 



(n = \4>\ 2 ), which is similar to the one used in the clas- 
sical O(N) model, no,fc denotes the condensate density 
in the equilibrium state. Note that we have introduced a 
second-order time derivative term. Although not present 
in the initial average effective action Ta, we shall see that 
this term plays a crucial role when d < 3 [461 US]. As 
pointed out in Sec. II D| the derivative expansion gives 
access only to the low-energy limit |p|,|u;|/c<Cfcof the 
correlation functions. It is however possible to extract 
the p dependence of the correlation functions by stop- 
ping the flow at k ~ (p 2 + lo 2 /c 2 ) 1/2 [8j. 

In a broken symmetry state with order parameter cf>i — 
s/2riQ, <p2 = 0, the two-point vertex is given by 



r fc 2 n(p) = v A,kU 2 + Z Atk e p + 2X k n Q ^ 
r fe 2 22(p) = ^A,fc^ 2 + Z A:k e p , 

,(2) 



(134) 



Using d82|, we then find 



S fc ,„(p) = G - 1 (p) + - 



1 r 



.(2) 



k :>n(p)+r k %( P ) 



• r (2) / \ ' - id 

u k,i2\P) fluid [13 . Equations (140 1 also imply that the Goldstone 



= A* + ^A,feW 2 + (1 - Z c , k )iuj 
- (1 - Z A , k )e p + A fe n ,fc 



(135) 



and 



(p) 



r!Cii(p)-rg a (p) =\ k n , k . (136) 



At the initial stage of the flow, Z a ,a = Zc.a = 1, Va.a = 
0, X\ = g and n 0j A — l^/g, which reproduces the results 
of the Bogoliubov approximation. 

Since the anomalous self-energy Sjt = o,an(p) ~ (^ 2 + 
c 2p2-)(3-d)/2 - g s i n g U i ar f or |p>| s joj]/c <C pg and d < 3, 
we expect ^ k>an (p = 0) ~ k 3 ~ d for k -C Pg (given the 
equivalence between k and (p 2 + uj 2 / c 2 ) 1 / 2 ), i.e. 



X k ~k 



3-d 



(137) 



The hypothesis (137) is sufficient, when combined to 



Galilean and gauge invariances, to obtain the exact in- 
frared behavior of the propagator. Furthermore we shall 



see that it is internally consistent. In the domain of va- 
lidity of the derivative expansion, IpJ 2 , \uj\ 2 /c 2 <C k 2 <C 
k 3 ~ d for k — > 0, one obtains from (134 1 



G k ,u(p) 

G k ,22{p) 

Gfe,i 2 (p) 



1 



1 1 

V A ,k + c 2 p 2 : 

Zc.k 



(138) 



2X k n 0tk V A M uj' 2 + c 2 k p 2 ' 



where 



Cfc = 



Z A , k /2m 



1/2 



V A ^ k + Z Ck /2X k no tk 



(139) 



is the velocity of the Goldstone mode. From (1371 and 



(1381, we recover the divergence of the longitudinal sus- 
ceptibility if we identify k with (p 2 + w 2 /c 2 ) 1 / 2 . 

The parameters Z A ^ k , Zc, k and V A<k can be related to 
thermodynamic quantities using Ward identities [5J H2J 



n s ,k = Z Atk n 0ik = n k 
1 d 2 U k 



Va 



k 



2n , fc dp 2 
d 2 U k 



(140) 



J C,k 



dndfi 



dn, 



o.fc 



dfi 



where h k is the mean boson density and n s ^ k the su- 
perfluid density. Here we consider the effective poten- 
tial U k as a function of the two independent variables 
n and fi. The first of equations ( 140 I states that in a 
Galilean invariant superfluid at zero temperature, the 
superfluid density is given by the full density of the 



mode velocity c k coincides with the macroscopic sound 
velocity 01131 IS], i.e. 



dn k 
dfi 



(141) 



Since thermodynamic quantities, including the conden- 
sate "compressibility" driQ ik /dfj, should remain finite in 
the k — > limit, we deduce from (140) that Zc, k ~ X k ~ 
k 3 ~ d vanishes in the infrared limit, and 



lim Cfc 

fc-S-0 



lim 

fc-fO 



( z A . k 

\2mV A u 



1/2 



(142) 



Both Z Aik = nfc/no,fc and the macroscopic sound velocity 
Cfc being finite at k = 0, V A , k (which vanishes in the 
Bogoliubov approximation) takes a non-zero value when 
k -> 0. 

The suppression of Zc,ki together with a finite value 
of V A , k =o shows that the effective action (133) exhibits a 



"relativistic" invariance in the infrared limit and therefore 



1G 



becomes equivalent to that of the classical 0(2) model in 
dimensions d + 1 |54| . In the ordered phase, the coupling 
constant of this model vanishes as Afc ~ ^-(d+i) ^ gee 
Sec. II D), which is nothing but our starting assumption 
(1371. For k — > 0, the existence of a linear spectrum 



is due to the relativistic form of the average effective 
action (rather than a non-zero value of Xk^o.k as m the 
Bogoliubov approximation). To neglect the term Zc,kd T 



in the average effective action ( 133 1 (and therefore obtain 
a relativistic symmetry) , it is necessary that Afc 3> k 2 [8] , 
a condition which is related to the singularity of the self- 
energies in the limit p — > 0. Thus we recover the fact 
that singular self-energies are crucial to obtain a linear 
spectrum in spite of the vanishing of the anomalous self- 
energy. 

To obtain the limit k = of the propagators (at 
fixed p), one should in principle stop the flow when 
k ~ (p 2 + ui 2 1 'c 2 ) 1 / 2 : . Since thermodynamic quantities 
are not expected to flow in the infrared limit, they can 
be approximated by their k = values. As for the lon- 
gitudinal correlation function, its value is obtained from 



the replacement 



C(w 2 +c 2 p 2 )( 3 - 



d)/2 



(with C a con- 
stant). From (1381 and (1401, we then deduce the exact 
infrared behavior of the normal and anomalous propaga- 
tors (at k = 0), 



name 



1 



n U)' 

r2 



c 2 p 2 



mc driQ icj 



n du lo 2 + c 2 p 2 



-G^ip), (143) 



Gan(p) 

where 



n^rac 



1 



n 



■ 2 ~ 7;Gn{p), 

uj a + c A p z 2 



Gn(p) = 



2n C(0J 2 + C 2p2)(3-d)/2 



(144) 



The hydrodynamic approach of Sec. IIIB correctly pre- 
dicts the leading terms of ( 143 ) but approximates dn^/d/i 
by n/mc 2 . On the other hand, it gives an explicit expres- 
sion of the coefficient C in the longitudinal correlation 



function (144 1 



2. RG flows 

The conclusions of the preceding section can be ob- 
tained more rigorously from the RG equation satisfied 
by the average effective action. The dimensionless vari- 
ables 



satisfy the RG equations 

3 ~ 1 - 
d t n ,k = - (d + r]c,k)no,k + ^Ik.u + -Ik,tt, 

<9 t Afc = (d - 2 + r) A ,k + ?7c,fe)Afe 

- A 2 , [9J fe;H) „(0) - 6J fc . tt>tt (0) + Jk:tt,tt( 

d - 

VA,k = 2X 2 k h , k —[j k . uu (p) + Jk;tt,ll(P) 

oy 

+ 2Jk;lt,lt(j>)] p=0 , 

g _ 

VC,k = - 2X k no t kg^[Jk;tt,lt(j') ~ Jk;tt,tt(p) 

- 3J k ;ii,it(p) + 3Jk;it,u(p)] p =o> 

d t VA,k = (2 - r\A.k + 2r]c,k)VA,k 


- 2\ 2 k h ,k-Q^[Jk-M,tt(p) + Jk-ttdi{p) 
+ 2J k ;itAt{p)] p= v 



(146) 



* 2 /fc 2 



where r} A ,k = -d t lnZ A , k , i) C ,k = -d t lnZ c ,k, V = P" 
and Cj = uiZc,k/ZA,k£k- The definition of the threshold 
functions / and J can be found in Ref. |5j. 

The flow of Afc, Zc,k an d V A ,k is shown in Fig. [6] 
for a two-dimensional system in the weak-coupling 
limit. We clearly see that the Bogoliubov approxima- 
tion breaks down at a characteristic momentum scale 
Pa ~ \J (gm) 3 n. In the Goldstone regime k <C Pd we 
find that both Afc and Zc, k vanish linearly with k in agree- 
ment with the conclusions of Sec. |III C 1| Furthermore, 
V A ,k takes a finite value in the limit k — > in agreement 



with the limiting value ( 142 ) of the Goldstone mode ve- 
locity. Figure [7] shows the behavior of the condensate 
density no,fe, the superfluid density n Stk — Za^Uq^ and 
the velocity Cfc. Since Z A , k =o — 1.004, the mean bo- 
son density h k = n s k is nearly equal to the condensate 
density no,fc- Apart from a slight variation at the be- 
ginning of the flow, no,fc, n Stk = Z Atk no tk and Cfc do not 
change with k. In particular, they are not sensitive to 
the Ginzburg scale pq. This result is quite remarkable 
for the Goldstone mode velocity c k , whose expression 
(139 1 involves the parameters Afc, Zc,k an d V A . k , which 
all strongly vary when k ~ pa ■ These findings are a nice 
illustration of the fact that the divergence of the longitu- 
dinal susceptibility does not affect local gauge invariant 
quantities [51 



3. Analytical results in the infrared limit 



ho t k = k Z c ,kno,k, 



7 -l 17-1 \ 



Afc — k' t k /. , ;./.,• A/, ■ 



(145) 



In the Goldstone regime k <C pa, the physics is dom- 
inated by the Goldstone (phase) mode and longitu dinal 
fluctuations can be ignored. If we take the regulator ( 132 1 
with r(Y) — ^T^O(y), the threshold functions I and J 
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\n(p G /k) 

FIG. 6. (Color online) X k , Zc,k and Va,h vs ^(pa/k) where 
pa = v (gm) s n/4iv for n = 0.01, 2mg = 0.1 and d = 2 
[ln(pG/Pc) — —5.87]. The inset shows vs 2m<? obtained 
from the criterion Va, pg = Va,/c=o/2 [the Green solid line is 
a fit to p G ~ (2mg) 3/2 ]. 

lr 




-5 

ln(pG/fc) 



FIG. 7. (Color online) Condensate density no,*,, superfluid 
density n 3i fc and Goldstone mode velocity c*, vs ln(pG /k). The 
parameters are the same as in Fig. [6] 
40 r 




ln(p G /k 



FIG. 8. (Color online) X' k vs ln(p G /fc) [Eq. |l50| ]. The pa- 
rameters are the same as in Fig. [6] 



can be computed exactly and one obtains [8] 

d t n 0tk = -(d + rjc,k)no,k, 

V d +1 X\ 



d t X k = (d-2 + T)c,k)X k + 8 



, v d+l Afe 

d+lv 1 ' 2 '' 

v A,k 



(147) 



dtV A ,k = (2 + 2ric,k)V A ,k, 



while r/A,k — 0. The first and last of these equations 
can be rewritten as riQ yk = no,fe=o and Vam — VA,k=o- 



From ( 147 I , we deduce 



d t \k = (1 - e)A A 



d t r]C,k 



-er)c,k ~ Vc,ki 



(148) 



where e = 3 — d. For d < 3, this yields Afe <~ fc^l — e) and 
lim 7]c, k = -e, (149) 

AC— 

i.e. AfejZc.fe ~ fc e in agreement with the numerical 
results of Sec. |IIIC 2| and the analysis of Sec. Ill C 1 



The anisotropy between time and space in the Goldstone 
regime k <S^pc (where the average effective action takes 
a relativistic form) can be eliminated by an appropri- 
ate rescaling of frequencies of fields. This leads to an 
isotropic relativistic model with dimensionless conden- 
sate density and coupling constant defined by [5] 



'0,/s 



VA,kno,k, 



A* 



(150) 



V A ,k 



X' k satisfies the RG equation 



d,Xi 



e\' k + 8 



d+l Ak ' 



(151) 



which is nothing but the RG equation of the coupling 
constant of the classical 0(2) model in dimensions d+1 
[Eq. (71 1], The corresponding fixed point value can be 
deduced from (74) [55J. In the infrared limit, we find 



X k = k- d (Z A . k e k ) 3/2 V^X' k ~ k% 



(152) 



if we approximate Z A , k ^ Z A>k=0 and V A ,k — Vk,fc=o- 
The vanishing of Afe ~ k e and the divergence of the longi- 
tudinal susceptibility is therefore the consequence of the 
existence of a fixed point X' k * for the coupling constant 
of the effective (d + 1) -dimensional 0(2) model which 
describes the Goldstone regime k <C pg- T° describe 
the entire hydrodynamic regime k <C p Cl we should in 
principle relax the assumption VA,k — fc=0j since Vajz 
strongly varies for k ~ po, which makes the analytical 
solution of the RG equations much more difficult. In 
Ref. |51| . it was shown that Eq. (151) is nevertheless in 



good agreement with the numerical solution of the flow 
equations in the entire hydrodynamic regime. We can 
then use ( 76 1 to obtain the Ginzburg momentum scale 



Pg 



.(d + l)e 



1A 



(153) 



in the weak-coupling limit, which agrees with the results 
of Sees. UTS] and MM 
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Appendix A: Threshold functions 



In conclusion, we have studied the classical linear 
0(N) model and zero-temperature interacting bosons us- 
ing a variety of techniques: perturbation theory, hydro- 
dynamic approach, large-iV limit and NPRG. We have 
shown that in the weak-coupling limit these two systems 
can be described along similar lines. They are character- 
ized by two momentum scales, the hydrodynamic scale 
(or healing scale for bosons) p c and the Ginzburg scale 
Pg- For momenta |p| <C p c , we can use a hydrodynamic 
description in terms of amplitude and direction of the 
vector field <p in the O(N) model, or density and phase 
in interacting boson systems. The hydrodynamic de- 
scription allows us to derive the order parameter correla- 
tion function without encountering infrared divergences. 
In the Goldstone regime |p| <C Pg, amplitude (density) 
fluctuations play no role any more and both the trans- 
verse and longitudinal correlation functions are fully de- 
termined by direction (phase) fluctuations. In this mo- 
mentum range, the coupling between transverse and lon- 
gitudinal fluctuations leads to a divergence of the lon- 
gitudinal susceptibility and singular self-energies. A di- 
rect computation of the order parameter correlation func- 
tion (without relying on the hydrodynamic description) 
is possible, but one then has to solve the problem of in- 
frared divergences which appear in perturbation theory 
when |p| < pg and signal the breakdown of the Gaussian 
approximation. The NPRG provides a natural frame- 
work for such a calculation. In the case of bosons, it 
shows that in the Goldstone regime |p|, \uj\/c <C pc-, the 
system is described by an effective action with relativis- 
tic invariance similar to that of the (d + l)-dimensional 
classical 0(2) model. 



The threshold functions appearing in the NPRG equa- 
tions for the 0(N) model (Sec. II D I are defined by 



I a = - fl(q)G£(q), 



(Al) 



4/3 (p) = <VJ a/3 (p). 



where a,/) e {l,t}. To alleviate the notations, we drop 
the k index. In dimcnsionlcss form, 

/>oo 

L = 2v d dyy d / 2 ( W + 2yr')G 2 a: 
Jo 

/>oo 

J Q(9 (0) = 2v d / dyy d ' 2 {r,r + 2yr')G 2 a G Pl 
Jo 



J«p(0) = 4^ / dy y d / 2 {[ V r + („ + 4)t/r' + 2y 2 r"]G 2 a 



d jo 

- 2(1 + r + yr')(r,r + 2yr')G 3 a }(1 + r + yr')G 2 , 

(A2) 



where 



Gi 
G t 



y(l +r) + 2\p 
1 



(A3) 



1/(1 + r)' 



and we have written the cutoff function as i?^(p) = 
ZkP 2 r{y) with y — p 2 /k 2 and r(y) a k independent 
function. For the theta cutoff function introduce in 
Sec. 



can 



II D 2 



= =—2-Q(l — y), and the threshold functions 



je computed analytically 



These strong similarities between the classical linear 
0(N) model and zero-temperature interacting bosons 
disappear in the strong-coupling limit. For the 0(N) 
model, this limit corresponds to the critical regime near 
the phase transition, which has no direct analog in zero- 
temperature interacting boson systems. The only ap- 
proach that one can hope to extend to strongly-correlated 
bosons is the NPRG. Recent progress in that direc- 
tion, based on the Bose-Hubbard model, is reported in 
Ref. m. 



where 



L = -8^ I I 



■Mo) 



J^(0)=4^I?, 



1 



1 + 2\p 



A 2 

d + 2! «' 



V 

d + 2 



A 2 A, 



A t = l. 



(A4) 



(A5) 
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